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Abstract: The Matrix String Theory, i.e. the two dimensional U(N) SYM with N = (8, 8) 
super symmetry, has classical BPS solutions that interpolate between an initial and a final 
string configuration via a bordered Riemann surface. The Matrix String Theory amplitudes 
around such a classical BPS background, in the strong Yang-Mills coupling, are therefore 
candidates to be interpreted in a stringy way as the transition amplitude between given 
initial and final string configurations. In this paper we calculate these amplitudes and 
show that the leading contribution is proportional to the factor gj x , where x is the Euler 
characteristic of the interpolating Riemann surface and g s is the string coupling. This is the 
factor one expects from perturbative string interaction theory. 
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1. Introduction 

The M = (8, 8) SYM on a cylindrical 2D space-time with gauge group U(N) (hereafter referred 
to as Matrix String Theory (MST)) is expected to represent in the strong coupling limit a 
theory of type II superstrings H, 0, [| . In the naive strong coupling limit the action reduces 
to the Green-Schwarz action for free closed strings of various lengths. More than that, it 
has been suggested that this theory describes a second quantized superstring theory Q (see 
also and the review article Q). In a couple of recent papers, ||, [l3|], it has been pointec 



out that the MST contains BPS instanton solutions which interpolate between different initial 



and final string configurations. In our previous paper, [13], we remarked that this could be the 
clue for a comparison with perturbative string interaction theory, not only from a qualitative, 
but also from a quantitative point of view. This is what we intend to develop in this paper. 

The main aim of this paper is to show that the MST in the strong coupling limit in 
the background of a given classical BPS instanton solution reduces to the Green-Schwarz 
superstring theory plus a decoupled Maxwell theory, and to compute amplitudes in such 
background. Since the latter interpolates between an initial and a final string configuration 
via a bordered Riemann surface X (which represents a branched covering of the base cylinder) , 
the amplitudes can be interpreted in a stringy way as the transition amplitudes between two 
such configurations. We show that their leading term is proportional to g7 x , where x = 
2—2h—b is the Euler characteristic of the Riemann surface of genus h with b boundaries, which 
characterizes the given classical solution. This is the result one expects from perturbative 
string interaction theory. 

The above derivation is contained in section |2| and is organized as follows. We first review 
the salient features of MST. Then we set out to compute the strong YM coupling limit g. 



i 



First we compute it for a classical BPS background (part of the elaboration is contained in 
Appendix [A]) . Then we expand the action about this classical configuration by splitting the 
quantum modes in two sets: the Cartan and the non-Cartan modes. It turns out that the 
non-Cartan modes can be neatly integrated out (Appendix [B|). What remains is a quadratic 
action for the diagonal (Cartan) modes. However these modes are not individually well 
defined fields on the cylinder. A field interpretation is possible if we lift them to the covering 
S. We show that, if we do so, we obtain the Green-Schwarz theory plus the free Maxwell 
theory on the world-sheet S. Afterwards, we pass to the calculation of the partition function 
of this theory and of the amplitudes mentioned above. In this regards a fundamental role is 
played by the Maxwell field zero modes. With a careful computation one can show that this 
amplitude is proportional to the factor g x = gj x announced above. 



2. The Matrix String Theory and its strong coupling limit 
2.1. Euclidean MST and the instanton background 



To start with let us summarize the results of [13|. MST is a theory defined on a cylinder C 
with coordinates a and r. Its Euclidean action is 

^2 




+ ite-D^e- + e+D w e+) + i a e T Ti{x\ e] , (2.1) 



where we use the notation 

1 1 

w = -(T + ia), w = -(r-ia), A w = A T -iA a , A € , = A T + iA a . 

Moreover X % with i = 1,... ,8 are hermitean N x N matrices and D w X l = d w X l + 
i[A w ,X 1 ]. 

Fww is the gauge curvature. Summation over the i,j indices is understood. 9 
represents 1 6 N x N matrices whose entries are 2D spinors . It can be written as 9 T = (9~ ,0^), 
where ± denotes the 2D chirality and 9~ , 6>+ are spinors in the 8 S and 8 C representations of 
SO (8), while T represents the 2D transposition. The matrices I\ are the 16 x 16 SO (8) 
gamma matrices. For definiteness we will write them as 

7i' 
,7i 

and 7i,7i are the same as in Appendix 5B of [15|. 



(2.2) 



The action (2T) has J\f = (8,8) supersymmetry. In we singled out classical super- 
symmetric configurations that preserve a (4,4) supersymmetry. In this configurations the 
fermions are zero, 6 = 0, and X 1 = for all i except two, for definiteness X % ^ for i = 1, 2. 
Introducing the complex notation X = X 1 + iX 2 , X = X 1 — iX 2 = X', the conditions to 
be satisfied for such BPS configurations are 1 

F w ^ + ig 2 [X,X} = (2.3) 
D w X = 0, D € ,X = 0. (2.4) 

1 Notice that one can obtain anti-instantonic configurations by choosing an opposite polarization for the 
(4, 4) broken supersymmetries. 



In jn| we found explicit solutions of these equations. Prom a mathematical point of view, 
( t^3j , 0) can be identified with a Hitchin system on a cylinder. Each solution of (|2.3|)(|2.4|) 
consists of two parts: a branched covering of the cylinder via the relative X characteristic 
polynomial and a 'dressing' factor. 

This can be seen by parametrizing the solution as X = Y~^MY and A w = —iY~ 1 d w Y, 
where Y takes values in the complex group SL(N, C) 2 and the matrix M will be explained 
below. The dressing factor is contained in Y while the branched covering is determined by 
M. 

Now, our purpose is to expand the action (2.1) around a generic classical BPS solution in 



inverse powers of the YM coupling g. Since, clearly, the background depends on g, we have 
to discuss preliminarily the strong coupling limit of the background itself. 

Let us review the branched covering part first. We consider the polynomial 

N-l 

P x (ji) = Det(> - X) = n N + M*Oi. 

i=0 

where fj, is a complex indeterminate. Due to (|2.4|), we have d w ai = which means that the 
set of functions {a^} are antianalytic on the cylinder. Therefore the equation 

Px(fi) =0 

identifies in the (w,fj,) space a Riemann surface S, which is an N-sheeted branched cover- 
ing of the cylinder. The explicit form of the covering map is given by the X eigenvalues 
set {x^\w), . . . ,x( N \w)}. The branch points are the locus where two or more eigenvalues 
coincide, which means where the identification cuts in the sheets start or end. As for the 
parametrization of the branched coverings we will choose the standard one 

/— ajv-i —a-N-2 — «o\ 

1 



M 



1 ... 



\ ...10/ 



(2.5) 



Notice that the branched covering structure is completely encoded in the {a{\ analytic func- 
tions and is independent of the value of the coupling. 

The dependence on the coupling is entirely contained in the factor Y . Unfortunately we 
do not have yet an explicit analytical treatment of the way the strong coupling configuration 
is reached by this factor, valid for any kind of covering. Therefore we limit ourselves to outline 
its features in general, extrapolating the validity of our discussion from the example of the 
%n coverings, which we are able to deal with more effectively, fL3|| . In Appendix [A| we work 
out some examples in detail. From those results it is sensible to assume that Y tends, in the 
strong coupling limit, to a precise matrix Y s (independent of g). This conclusion needs a more 
precise statement: if we cut out a neighborhood of size proportional to some positive power 
of 1/g around each branch point, then YY^T 1 = Yd dies off to 1, outside such a neighborhood, 
more quickly than any inverse powers of g. Yd will be called the dressing factor. Since in this 
2 See Appendix |a] for the change from GL(N, C) to SL(N, C) with respect to Q 



paper we are interested in expanding the action Q2.1|) in inverse powers of 1 /g, and actually 
in singling out the dominant term in this expansion (see below), we will consider the action 
( |2,1| ) around a given classical solution stripped of the above dressing factor and exclude from 
the integration region the branch points on the cylinder. In other words we will consider from 
now on the action ( j2.1|) in which the pertinent Y is replaced by Y s and the integral extends 
over Co which is the initial cylinder C from which the branch points have been removed. In 
the Ztv case the surviving factor Y s is (V J J) -1 , with the notation of ||13|| . 

After getting rid of the dressing factor, the classical background configuration is specified 
by X = Y~ 1 MY S and A w = —iY~ l d w Y s . As expected, this configuration is singular exactly 



at the branch points. Now one can diagonalize the matrix M (2.5), M = SXS , with 
a nonsingular matrix S which can always be decomposed as the product of a Hermitean 
matrix W and a unitary matrix T, S = WT, [T^l , (in the Zjv covering case, W = J -1 and 
T = A -1 , see fl~3|l ). X is the matrix of eigenvalues of X and X = Y S ~ 1 SXS~ 1 Y S in the strong 
coupling limit. Remarkably, Y~ 1 S is itself a (singular) unitary matrix and so simultaneously 
diagonalizes X and X. Corresponding to X we have A w = —iS~ l d w S which can be seen 
to be zero everywhere, even at the branch points. What has happened is that the unitary 
transformation has swallowed entirely the connection, including the singularities. 

We have therefore two distinguished (but equivalent) ways to represent the classical back- 
ground in the strong coupling limit: 



• strong coupling representation: X s = Y~ 1 MY S and A s = —iY~ 1 dY i 

• diagonal representation: X is diagonal and A = 0. 

The first is the natural form of the background on the cylinder once the dressing factor is 
removed by the strong coupling limit. The second, as we will see later on, is suitable to 
represent the theory on the covering space E. 

To end this subsection we make two comments. The first concerns the consistency of our 
procedure. One may feel uneasy for the appearance of singularities and the use of singular 



gauge transformations, as above. As stressed in [13], the classical configuration specified 
by a given couple (X, A) is smooth on the initial cylinder C, but if we strip the solution 
of the dressing factor we get a configuration which is singular exactly at the branch points. 
The dressing factor is there exactly to compensate for these singularities. As pointed out in 
Appendix [A], Yd — 1 has support only at the branch points in the strong coupling limit. It 
is perhaps possible, but formally very complicated, to keep such factor in the action. In this 
paper we prefer to replace C with Co by excluding the branch points from the integration 
region to preserve smoothness. This is justified by the following consideration. Beside the 
initial smooth configuration (X, A) on C, we will meet another smooth situation when we lift 
our theory to the branched covering £ of C (see below). At that stage the branch points can 
safely be restored in the integration region. What has happened is that, in order to pass to 
the covering, we need a singular gauge transformation, the Xs -1 !? used above, which exactly 
kills the singularity exposed by the strong coupling limit. In other words, it is natural to 
perform this singular gauge transformation if we want to reach a smooth situation which is 
fit for field theory. 

We can regard the same problem from a different perspective. The moduli space of couples 
(X, A) satisfying (OIO) have b een studied with sophisticated mathematical methods ]l0|. 
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However we can say, roughly speaking, that they consist of the group degrees of freedom (the 
factor Y and their generalizations, modulo gauge transformations) times the moduli space 
of the Riemann surfaces determined by the branched coverings. The strong coupling limit 
suppresses the former and only the latter survive. Since the risk of using singular gauge 
transformations is to suppress degrees of freedom or to introduce new ones, we see that in 
our case this does not happen: the moduli present in the theory after lifting it to the covering 
will correspond to those that have not been suppressed by the strong coupling limit. 

The second comment concerns another way to get the strong coupling limit of the back- 
ground; here we simply outline it. In the parametrization 

A w = -iY^dyjY , X = Y~ X MY , 

we get 



Fww — d w A w d w A w + i [A w , A w ] — %Y d w d w (YY^^ (YY^^ 
and 

[X,X] = Y- 1 \m, (YY + ) M f (YY+y 1 ] Y . 



Y 



Therefore the BPS equation on Y is 



d w (d^nn- 1 ) + g 2 



o, 



where $7 = YY + . Notice that this equation is the extremality condition for the deformed 
WZNW action 

i \n] = — / d 2 wTvn~ 1 d w nn~ 1 d w n+ 

167T J 

X — J d 3 ye ijk Tr (jr^nfi-^nfi -1 ^) - |- J d 2 wTr (jr 1 MftM + ) . 



24vr JB\dB=Cyl 

One could follow the g deformation of the space of classical solutions of the above theory 
subject to appropriate boundary conditions and obtain its large g limit analytically. 

2.2. Fixing the gauge and integrating along the non Cartan directions 

To extract the strong YM coupling effective theory, we first rewrite the action in the following 
useful form 

S = ^Jd 2 wTr (d^D^X 1 -^-[X i ,X j ] 2 -g 2 [X I ,X][X I ,X]+D w XD w X 

-^2 (F ww + ig 2 [X,X}) 2 + ifoDvOj + 9+D w 9+) + igO^iX 1 , 0]j , 

where / = 3,4, ...,8. We now expand the action around a generic instanton configuration 
writing any field <5 as 

$ = $( fc ) + ^ + n = $( 6 ) + ( /, = $ o + 0", (2.6) 

where is the background value of the field at infinite coupling, (p 1 are the fluctuations 
along the Cartan directions and <p n are the fluctuations along the complementary directions 
in Lie algebra u(N). Of course only the upper case fields X, X and A will have non zero 



background value. By their background value we mean either the strong coupling or the 
diagonal representation introduced above. In the latter case (j) 1 and eft" are the diagonal and 
non-diagonal part of (f>, respectively. In the strong coupling representation cfi 1 is the part of <f> 
such that Tr(c/>X™) 7^ for some n, while (fi n is the part of <p for which Tr(</>A") = for any 
power n > 0. Anything we say in this subsection holds for both representations. 

The first important remark is that, since the value of the action on the background is zero 
and since the background is a solution of the equation of motion, the expansion of the action 
starts with quadratic terms in the fluctuations. To proceed further we have to fix the gauge. 
We use, in the strong coupling limit, a gauge fixing inspired by t'Hooft 'abelian' gauges, and 
similar to the one used in [14]. We write it in the form 

G wiD = D° w a w + D%a w + ig\[X°,x) + [X°,x]) + 2ig 2 [X oI ,x I ] =0, (2.7) 

where D° stands for the covariant derivative with respect to A°. Next we introduce the 
Faddeev-Popov ghost and antighost fields c and c and expand them like all the other fields. 
Then we add to the action the gauge fixing term 



S 9f = 4^2 / ^Iw « 2 '' 



and the corresponding Faddeev-Popov ghost term 



1 f ,2 -&Gww 



S ghost = "2^2 J dW ( 2 ' 9 ) 

where 5 represents the gauge transformation with parameter c. 

At this point, to single out the strong coupling limit of the action, we rescale the fields in 
appropriate manner. Precisely, we redefine our fields as follows 

A W = A$ +gal + al, X = X® +x l +-x n , X 1 = x 11 + -x In , 6 n = 9 l + ^-9 n 

9 9 y/9 

and likewise for the conjugate variables. For the ghosts we set 

c = gc l + -y/gc n , c = gc t +^-c n . 

These rescalings introduce a unit Jacobian in the path integral measure of the non-zero 
modes, but they may produce a non-trivial factor due to the presence of zero modes (see 
below). 

After these rescalings the action becomes 



S = S sc + Q n + o 

where 



1 
71 



S sc = - [ d 2 w Tr \D^x n D^x n + D^x^h' + i{elD^el + elD^el) 

* J Co l 



+D^a t -D { - ) a t + D^c^c 1 
~ru w u> w u w u, w ~r -is,,, <*■ J-^ii) <-- 



(2.10) 



and represents the covariant derivative with respect to the background connection. Q n 
is the purely quadratic term in the (j) n fluctuations. The latter can be easily integrated over 
and, since they do not involve zero modes, give exactly 1, in accord with supersymmetry (see 
Appendix |l|) . 

In conclusion, in the strong coupling limit we are left with the quadratic action (|2.1C| ) 
over the Cart an modes. 



2.3. Lifting to the branched covering 



Let us now show that the effective theory we obtained in the previous subsection corresponds 
to the Green-Schwarz superstring plus a free Maxwell action on the worldsheet identified by 
the branched covering of the relevant background. To this end we consider the quadratic 
action ( |2.10| ) in the diagonal background: the covariant derivative D$ becomes the simple 
derivative and the action becomes a free action 



s S c = - [ d 2 w Tr id^d^x 11 + d^d^x 1 + iield^l + e^el) 

* J Co 

+d w a t iB du,a i w + d^d^c 1 ] . (2.11) 



Since all the matrices are diagonal we can rewrite this action in terms of the diagonal modes 



^ = ~ j c d 2 W Y^n=l 9 v> x \n) d ™ X \n) + ] s(n)dw6 s(n) + dc(n)d w 9 c (n)) 

(2.12) 



This is a theory of free fields on Co and it is tempting to extend the action to C by just 
forgetting the punctures on the cylinder corresponding to the branch points. However this 
is not correct. The fields x % are not single-valued on the cylinder. For example, in the Zjy 
covering case, upon going around a branched point, each x % is mapped to the adjacent one, 



and this is precisely the reason why they describe long strings, ||[13]. Mathematically, this 
problem can be rephrased as follows: the fields in ( |2.12| ) are not section of bundles over C but 
they can be suitably combined to form sections of line bundles on the covering £ of C. 

At this point it is worth spending a few words about Hitchin systems. The Hitchin systems 
we are interested in are defined starting from a U(N) vector bundle V over C, associated with 
the fundamental representation of U(N). They consist of couples (A, X) where A is a gauge 
connection and X a section of EndV (8> K , where K is the canonical bundle of C, which satisfy 



(fO|) and Q2.4D, [1C]. Such systems can be lifted to an A r -branched covering of C , [ 1 1 1 , [ 16 j , IT 



A remarkable feature of the lifting is the appearance on the branched covering of a line bundle 
L constructed out of V and from which in turn V can be reconstructed. In simple words 
the initial non-Abelian system can be described by an equivalent Abelian system on the 
branched covering. This is exactly the situation we are faced with when lifting our action to 
the branched covering. Without using rigorous mathematics, let us try to render this fact 
plausible by looking at the realization of local fields on a Riemann surface represented as a 
branched covering. If S is a branched covering of the cylinder C, then there exists a projection 
map 7r : S — > C whose inverse image is N-valued. In our language this is simply 

vr -1 : w -> (i( 1 )H...,i(jv)(tt)). (2.13) 
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Suppose a local complex field tp is given on E; applying the above construction if) can be 
represented as an N-tuple (jjjru(w), . . . ,i/j(N)i w )) representing the field on each copy of the 
cylinder C that composes the covering E; the ^)(ty)'s are related by the appropriate mon- 
odromy properties along the curves $tw = const. 

Going now back to the action (|2.12j ), we will interpret any set of N fields (4>(i)-, ■ ■ ■ , <^(jv)) 
in it as a unique field <j) on the covering E. (p is locally a function of a coordinate z in E. From 
the point of view of E, the w coordinate is locally defined via an abelian differential to = dw 
with imaginary periods which is canonical, i.e. is fixed only by the complex structure of the 
surface (see @ and 01). 

Finally we can write the strong coupling action ( 2.12j ) as follows 



Sfc — + ^Maxwell' (2-14) 

s gs = IJ s d 2 z (d z z l d- z x l + i(6 s d z e s + e c d z e c )) (2.15) 

- f d 2 z(g zS d z a- z d- z a z + d z ~cd- z c) . (2.16) 
71" Js 



'-'Maxwell 

7T 



In ( |2.15| ) a ^fuJl (resp. t/ujz) factor has been absorbed in 9 S (resp. 9 C ) which is a (5,0) (resp. 
(0, g)) differential on E and the metric in the Maxwell term is g zz = oo z LOz- 

Summarizing, what we obtained in this subsection is that the strong coupling effective 
theory is given by the Green-Schwarz superstring action on the branched covering worldsheet 
plus a decoupled Maxwell theory on the same surface. 

The fields in ( gig ) are now sections of line bundles over E, i.e. well defined fields on the 
Riemann surface: for example a is a section of the canonical bundle of S and so on. 

2.4. String amplitudes 

Let us compute first, for simplicity, the vacuum to vacuum amplitude of the SYM theory 
in the strong coupling limit in the background of a given instanton. As we have already 
pointed out several times, this amplitude (up to the vertex string insertions, see below) has 
a string interpretation as the amplitude for the transition from the initial to the final string 
configuration described by the instanton. If this interpretation is correct this amplitude, to 
the leading order, should be proportional to gj x where x 1S the Euler characteristic of the 
Riemann surface E, i.e. the covering surface introduced above. 

What remains for us to do in order to evaluate this amplitude is to integrate over the 
Cartan modes in the functional integral with action ( p,15|) (the non-Cartan modes have 
been integrated out above). Since the action is free, the integration produces a ratio of 
determinants, which turns out to be a constant. However we have to take account of the zero 
modes for the fields that have been rescaled (the unrescaled zero modes are irrelevant in this 
argument). The rescaled fields in C are the Maxwell and the ghost fields. The corresponding 
fields in E will be rescaled too 

a z ^ga z , a z ^ga z , c^gc, l^gl. (2.17) 

Therefore let us single out the Maxwell (plus ghost) partition function. We will show that 
the decoupled U(l) theory is there to generate the stringy factor gj x as a consequence of the 
rescaling (2.17), (on the Maxwell partition function, see also section H). 



In fact under this rescaling the Maxwell partition function (a = a z ,a = a 2 ) 
Zl axweU = Jv[a,E,c,~c] e -*S« we «(^,a,g) 

rescales with a factor depending on the zero modes. Roughly speaking, what happens is that 
the above integral is interpreted as the ratio (Det / D c )/(Det'n a ), where □ = dd denotes the 
quadratic operator in the action, and ' means that the zero modes have been excluded from 
the computation of the regularized determinants; since we have rescaled the measure there 
will arise a factor of g at a power equal to the unbalance of the zero modes 3 . The problem 
is therefore to count the latter. As for the ghost fields which are scalars, the only zero modes 
of the d operator on E is the constant. The zero modes of the Maxwell fields correspond 
to the holomorphic differential on S. If E were a closed Riemann surface of genus h, their 
number would be h. Their counting in the present case is not completely standard as E is 
actually a Riemann surface with boundaries (representing the in- and out- strings). A way 
to do the counting is to construct the double S of S: E has genus h = 2h + b — 1, where 
b is the number of boundaries, and admits an anticonformal involution with the set of fixed 
points corresponding exactly to the boundary of E. We can count now the number of analytic 
differential on E that extend to S, that is the so-called analytic Schottky differentials |p0| : 
their number is h. Therefore the overall unbalance of zero modes (including the ghosts) is 
h — l = 2/i + 6 — 2, which is exactly the opposite of the Euler number of E. An equivalent 
way of deriving this result is to use the Gauss-Bonnet theorem on E and noting that, due to 
the involution, the integral of the curvature over E is one half of the total contribution. 

-2/i-b+2 _ n 2h+b-2 



Finally the factor in front of the vacuum to vacuum amplitude will be g + = g 



.s 



The exponent of g is precisely the Euler characteristic of E, as we wanted to prove. 

In order to appreciate exactly what we have just computed we must now specify what it 
corresponds to in string interaction theory. In this sense the amplitude we have just computed 
in the strong coupling limit is a basic amplitude but, of course, an incomplete one. 

First of all, real string amplitudes should contain vertex string insertions, i.e. should be 
correlators of the vertex operators corresponding to the various in- and out- (super)strings. In 
this regard we simply remark that such vertex operators are constructed in terms of the string 
fields x % and 9 (and, possibly, the string ghosts and superghosts), therefore the treatment of 
the non-Cartan modes above is not affected and the discussion of the zero modes of the 
Maxwell sector is unchanged. Therefore the scaling factor gj x is left unchanged too. 

Moreover, in order to obtain complete amplitudes, we must still integrate over the moduli 
of the Hitchin systems, i.e. over the inequivalent BPS instantons that interpolate between a 
given initial and a given final state. If we want to implement such more advanced stage of cal- 
culation, we have to take into account in the measure some Jacobians that are produced by the 



various field splittings we have considered above. Following [22], the background/fluctuations 
splitting of the fields in the path integral generates a Jacobian Jb/f- In an analogous way also 
the Cartan/non-Cartan splitting gives rise to another Jacobian factor Jc/nC- These factors 
are easily seen to depend only on the Cartan modes of X and 6. This, in particular, implies 
the validity of the procedure we used for the integration over the non-Cartan modes. 



3 A more precise account of this point can be found in [|19j where an analogous calculation of the Maxwell 
partition function was carried out. 



n 



The total partition function in the strong coupling limit is therefore 

Z sc = I dm gjx I V[x, §, a, c] J b/f J c /nC^ SGS ~ SMaxw M , (2-18) 

J Msc J 

where Ai sc is the strong coupling instanton moduli space. But we know what this is: the 
strong coupling limit of each BPS configuration is completely identified by the relevant 
branched covering. This means that the integral over the background instantonic configura- 
tions corresponds in the strong coupling to the sum over all inequivalent Riemann surfaces 
with appropriate boundaries. What is important to stress here is the g^ x factor in front of 
the path integral at fixed genus. 

At this point the situation is clear: Matrix String Theory in the strong coupling limit looks 
very much like IIA theory. We would like to say 'is' instead of 'looks like'; however a complete 
proof of the equivalence of the two theories requires showing that the above Jacobian factors 
give rise to the full superstring measure in the path integrals over the surfaces |23[j ; this in 
turn requires a manageable representation of the moduli space of the Hitchin systems in the 
strong coupling limit. 



3. Comments 



In this final section we would like to make a few comments, in part to complete the previous 
presentation and in part to outline future developments. 

The first comment concerns the explicit breaking of the SO(8) invariance of ( |2.1[) in the 
background: we recall that, at the moment of choosing a fixed BPS background, we set X 1 7^ 
only for i = 1,2. As one can see, in the strong coupling limit no remnant of this explicit 
breaking survives. However the problem of SO(8) symmetry breaking may arise in non-zero 
orders of the 1/g expansion. In such a case SO (8) invariance is recovered by summing over 
the instantons lying along all couples of the eight directions. 

Next let us briefly discuss the relation of the results of this paper with the large N limit, 



24]. It must be stressed that what we have found above is in agreement with the present 
understanding of the M(atrix) theory pq] . In fact in order to see string worldsheets, and 
make computations with them, we do not need to take the large N limit: they are clearly 
visible and effective for finite N and the large ./V limit does not add anything to such visibility. 
Of course the U(N) theory describes only a limited amount of string processes due to the fact 
that the total content of each asymptotic sector, given by a definite long string configuration 
(l) ni • (2) n2 • . . . • (s) ns , has finite size N = Xli=i^ n «- The full infinite set of possible string 
processes can be truly obtained only in the large N limit. 

We would like to end the paper with some comments concerning the free Maxwell action 
and the related partition function, which appears in the strong coupling limit. Above Maxwell 
field a on S appears as a 'small' fluctuation. We would like to be more precise about this 
point. 

Following 1 29], the quantization of the Maxwell sector is straightforward, 

n 

4»e„W = E ( dim II x{RUe-^ {R) , (3.1) 



1 n 



where m = (mi, . . . ,m n ) are the periods of the U(l) connection around the boundaries of E, R 
is any irreducible U(l) representation, x(R)m 18 the relative character, = J E ^fgd?z is the 
area of E in the branched covering metric and C2OR) is the Casimir of R. 
Explicitly evaluating fl3.1|) , we get 



Z MaxwellW = 2^ e '' e 2 

lew 

(3.2) 

where |m| = a = J" 2 / and f = da is the curvature of a = a^dz + d^dz. We assume 
|m| to vanish, and this is our definition of Maxwell fluctuation in this paper. This means in 
particular that the Maxwell field is just a section of the canonical line bundle, as we assumed 
implicitly above while doing zero-modes counting. With this assumption the Maxwell factor 
is the constant 

S ^ NA Cyl ; 2 

Z Maxwell ~ e 2 ' 

lew 

which goes to 1 in the large N limit. 

We would like however to point out that the Maxwell theory on E admits other sectors 
in which |m| ^ 0. These non-trivial sectors have not been considered in the present paper. 
However their inclusion in the theory looks extremely interesting. We recall that, in a non- 
interacting regime (trivial background) of MST, the number of units of quantized electric flux 
is identified with the DO-brane charge, |26| , p7f| . The number |m| may have to do with the 
conserved DO-brane charge involved in the interaction. Related to this is the remark that, 
although the Maxwell sector is completely decoupled from the string sector of the theory in 
the strong coupling limit, 0(1/ g) and higher corrections to the effective action will bring into 
the game interaction terms representing non-perturbative contributions to string theory. 

A. Solutions for the Toda dressing factors 

In this Appendix we discuss the properties of the Toda dressing factor e u ' H in the classical 
instanton solutions in the strong coupling limit. We recall that this factor appears in the 
dressing matrix Y. With respect to [13] we introduce here a convenient modification: instead 



of matrices Y, J, U in GL(N, C) we will consider matrices Y, J, U in SL(N, C); simultaneously 



the matrices H and TL will be traceless. Any result can be obtained directly from [13] by 
simply dividing by the corresponding determinants or subtracting the corresponding traces. 
The reason for this is that the determinant/trace factors introduce a useless complication for 
our analysis. With this understanding we will use the same symbols as in |l3|| : in particular 
(Hi) kl = (5 hk - 6i, k +i)S k i, i = 1, ...N - 1 and k, I = 1, ...,N, and H = Diag((l - N)/2N, (3 - 
N)/2N,...,(N - 1)/2JV). 

In the case of the Z^v covering X N — a = the 'fields' u ■ H satisfy the Toda-like equations 

d ? 3 c -(u ■H)+g 2 \e- u - H Pe u - H ' ,P*\ = irH6(a)(d c a)(d£a) , . (A.l) 
where j$- = a l / N . In particular, in the case = 2, we have only one field u satisfying 



7T 



d c d^u - 2 5 2 sinh u = --5(a) (d c a) (d^a) . (A.2) 



1 1 



This means that u must satisfy the sinh-Gordon equation with the boundary condition 

u~— -ln|a|, a~0. (A. 3) 

Let us suppose, for simplicity, that a = z — zq, where zq is a point away from the origin of 
the z-plane, and z = e w . This corresponds to the case of N closed string joining to form one 



long string, see |13| . Now we want to single out, for this simple case, solutions of ( A. 2 ) with 
the right asymptotic behaviour to interpolate between the in- and out- string states. First 
we rescale ( — > V%gC, so that the sinh-Gordon equation takes the standard form 

d^d^u — sinh u = . (A. 4) 

We do not know an exact solution of this equation satisfying the boundary condition ( |A.3j ) 
and vanishing at the origin and at infinity of the z-plane. We can therefore proceed in two 
ways: find a numerical solution or an approximate analytical one. Let us do the latter first. 
Recalling that 



d 7 f = V2 9 ^^, 

oz z 

the approximate expression of ( in terms of z is 

C ~ '~7r^-{z ~ z of /2 , for z ~ z 
3z 

C r ^-V^ for \z\ » \zq\ 

C ~ V2g^/— zqIuz, for \z\ « \zo\ ■ 

If these were the exact expressions for we could consider spherically symmetric solutions 
of ( |A.4[ ), i.e. solutions depending only on r = For them eq.( |A.4[ ) takes the form 

d^u + -d r u = 4sinhn (A. 5) 

r 

This is a form of the Painleve III equation. The general form of the solutions of this equation 



are known, see [28] and references therein. Let us select the class of solutions with the 



following asymptotic behaviour: 

u(r) ~ alnr + (3, r -> 0, |a| < 2 (A.6) 
u(r) ~ 'yr~ 1 / 2 e~ 2r , r — > oo . (A. 7) 

The constants (3 and 7 must be fine-tuned to a in order to give rise to smooth solutions. 
However, in our case, we are not interested in the actual value of (3 and 7, therefore we can 
always adjust the parameters in such a way as to have a smooth solution. The only serious 
problem might come from the bound \a\ < 2; however the asymptotic behaviour ( |A.6| ) is in 
our case 

u ~ lnr 

6 

therefore the bound is satisfied. 



1 o 



Figure 1: Numerical solutions of eq. ( A.2) for (a) small and (6) large g, respectively. For the 
Bullough-Dodd case, eq. (A. 10), the solutions are almost identical. 



Let us study now the properties of the solution. In the various regions we have the 
following asymptotic expressions: 

2V2g . 



3|^o| 



zq| 3 ^ 2 , for z ~ zq 



r ^ _^£y^[ ? for \ z \ >> \ ZQ \ (A. 



V^^a/I^oI 111 !- 2 !) for ,\z\ « \zq\ . (A. 9) 



From this we see that, when, at fixed finite g, z is near the origin and far to infinity, the 
solution tends to zero. The convergence to zero is more rapid the larger g is, the slope in g 



being of negative exponential type. Looking now at the first equation ( AJS), we see that, even 



if we sit near zq, we may still fall in the regime (r large) in which the solution is extremely 
small, provided g is large enough. In other words, for large g the solution shrinks around zq, 
and, in the g — ► oo limit it becomes spike-like with support at z = Zq. We can say that, if we 
exclude a neighborhood of zq of size proportional to (g)~ 2//3 , the solution decreases to zero 
more rapidly than any power of 1/g. 

We recall that the spherically symmetric solution is not the exact solution, but only 
an approximate one. However we expect the general behaviour of the true solution to be 
essentially similar, i.e. that it shrinks very rapidly around zq as g —* oo. This is confirmed 
by the numerical solutions shown in fig. [|. 

We can easily extend the previous analysis to the case in which a contains several distinct 
zeroes. Simply find suitable approximate expressions for Q near the zeroes of a and apply the 
previous approximate analysis. The conclusion will be that the solution shrinks very quickly 
around the branch points as g — » oo. 

Let us briefly consider the case N = 3. We have two fields u\,U2 obeying the equations 

d ( d^ui + g 2 ( e - Ul ~ U2 - e 2ui ' U2 ) = ~5{a) d^a c\a 

d c d- c u 2 + g 2 (e- Ul ~ U2 - e 2u2 ~ Ul ) = --5(a) d^a d ( a . 

3 

Subtracting one from the other we get an equation which is identically satisfied with u\ = 
The remaining equation for u = ui = u 2 is 

d^u + g 2 (e~ 2u - e u ) = --5(a) c\-a d ( a . (A.10) 



1 -t 



This is known as the Bullough-Dodd equation. The features of the solutions, with the present 
boundary conditions, are essentially the same as for the sinh-Gordon equation: for g — > oo 
they shrink around the zeroes of a. The numerical solutions do not differ significantly from 
the sinh-Gordon case. 

We will not discuss the N > 3 case here. 



B. Non— Cartan fluctuations 

In this Appendix we briefly discuss the Q n term of the strong coupling action, i.e. the 
quadratic terms in the non-Cartan modes. Q n has the form 



Qn 



1 



7T 



where 



x n Qx n + x In Qx In + a%Qal + c n Qc n + i(6^, 9")A 



9" 



(B.l) 



and 



Q = adjf o • adx° + ad a t • ad a * + &d x it • ad^i 



A 



i&d a t_ 7;ad X c 
7iadjf ol iad a t 



In the path integral we can now integrate over the non-Cartan modes and obtain a ratio of 
determinants of A and Q. Since these operators do not have zero modes the calculation is 
elementary. The integration over a n and the conjugates exactly cancels the integration over c" 
and the conjugates. What remains is a ratio ((Det.A) 16 /(DetQ) 8 ) . The expression of the 
numerator is formal: one should understand Det.A as y / Det(— AA?). But AA^ = A^A = — Q. 
Therefore the net result of integrating over the non-Cartan modes is 1. This is the result 
expected from supersymmetry in the absence of zero modes. 
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